Optical sum rules that relate to the potential energy of strongly correlated systems 
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A class of sum rules for inelastic light scattering is developed. We show that the first moment 
of the non-resonant response provides information about the potential energy in strongly correlated 
systems. The polarization dependence of the sum rules provide information about the electronic 
excitations in different regions of the Brillouin zone. We determine the sum rule for the Falicov- 
Kimball model, which possesses a metal-insulator transition, and compare our results to the light 
scattering experiments in SmB6. 
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One of the key signatures of strong electron correla- 
tions in condensed matter systems is the redistribution 
of spectral weight from low to high energies. A clas- 
sic example of this phenomenon occurs when a normal 
metal becomes a superconductor at low temperatures. 
The optical conductivity is suppressed at low energy due 
to the presence of the superconducting gap Q . Because 
there is an optical sum rule — the /-sum rule — which con- 
strains the integrated spectral weight in cr(i>), the spec- 
tral weight suppressed below the superconducting energy 
gap reappears as a delta function peak at w = 0, reflect- 
ing the onset of many-body coherence in the system. In 
high temperature superconductors, low-energy spectral 
weight has also been observed to shift to high energies 
below the superconducting transition Q . This dramatic 
spectral weight transfer is believed to be a signature of 
the strong electronic correlations in these and related ma- 
terials 0. 

The use of sum rules in optical conductivity measure- 
ments has had a wide impact in a number of fields of sci- 
ence. The optical sum rule originated in atomic physics 
as a relationship between the total number of electrons in 
the atomic system and the integrated spectral weight . 
In a solid state system, the optical sum rule is usually 
projected onto the lowest energy band. In this case, the 
integrated spectral weight associated with electrons in 
the lowest band is related to the average kinetic energy of 
the electrons, regardless of the nature of the electronic in- 
teractions [f| . Since the electronic kinetic energy usually 
varies on an energy scale on the order of electron volts, 
the average kinetic energy is essentially a constant at low 
to moderate temperatures, and so the optical conductiv- 
ity sum rule is essentially temperature-independent. 

It would be useful to have similar sum rules that re- 
late to the potential energy. In strongly correlated elec- 
tron systems in which electronic interaction energies are 
on the order of, or greater than, the kinetic energy, ef- 



forts to study the potential energy as functions of dop- 
ing, pressure, or temperature have been hindered by the 
limited tools available for measuring dynamical proper- 
ties. We show here that the non-resonant inelastic light 
scattering response has sum rules similar to the f-sum 
rule that relate to the potential energy, and that these 
sum rules should provide valuable insight into the scat- 
tering of optical or x-ray photons from materials 0. As 
an added benefit compared to optical conductivity mea- 
surements, inelastic light scattering allows one to probe 
the electronic excitations in different regions of the Bril- 
louin zone (BZ) by orienting the incident and scattered 
light polarizations. 

Shastry and Shraiman [?J were the first to point out 
that, if the self energy is local and d — > oo, the optical 
conductivity and the nonresonant Raman response func- 
tion in the depolarized (ej perpendicular to eg) scattering 
geometry are related according to a(v) cx S{y)/v a rela- 
tionship that was proved in Ref. 0. This result implies 
that there is also a sum rule for Raman scattering, which 
is related to the f-sum rule for optical conductivity and 
to the average kinetic energy of the electrons. 

In this paper, we show that there are also interaction- 
dependent sum rules for inelastic light scattering that are 
related to the potential energy for q = Raman scat- 
tering, and to combinations of the potential and kinetic 
energies for q ^ inelastic x-ray scattering. 

The formalism used to derive these sum rules is 
straightforward. If we consider the time-ordered prod- 
uct that yields the susceptibility corresponding to the 
fluctuations of an operator O 

X (r) = TT(T T e-^ H O(r)OH0))/Z - (Tr(e^ n O)/Z) 2 , 

(1) 

where r is the imaginary time, 7i is the Hamiltonian, (3 
is the inverse temperature, Z = Tre~@ n is the partition 
function and 0(t) = e T7i Oe~ TH . Then it is easy to show 
from the spectral formula that when the susceptibility is 
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Fourier transformed, and analytically continued to the 
real frequency axis, we get 



dwJm X (q, v) = -Tr(e-^[C, [H, 0^}})/Z (2) 



as the general sum rule for the first moment of the suscep- 
tibility (the brackets indicate commutators). If we con- 
sider the polarization operator O po i = e J^. Rin^ with n, 
the number density operator at site i, then this gives 
the /-sum rule where <j(v) — \v&xjj{y)j'hv in terms 
of the current-current correlation function XJJ (with 
J = i[H,O pol ]/h: 
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-T?(e- m {-f))/Z. 



(3) 



The sum rule is proportional to the average kinetic 
energy (we consider only nearest neighbor hopping on 
a hypercubic lattice T = —t ^\ s <j c \ a- c i+s,cr with 5 a 
nearest- neighbor translation vector). Likewise, if we 
consider the dynamical charge susceptibility c h(q) = 
Efe c k+k/2,<T C k- q /2, ( r we obtain 



/ dwlmxch{q_, v) 
Jo 
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■-Tr(e-^J2i e ^- 



q) + e(k-q)-2e(k)])/Z, 



= 7r(l-X)Tr(e~ pAt (-T))/,2 ) (in d -> oo) (4) 

with X = linid^oo X^=i cos(ki)/d and e(q) the band 
structure. This vanishes as q 2 for q=0 as it must, due to 
total charge conservation. 

There is, however, no charge conservation for inelas- 
tic light scattering with q = optical photons (Ra- 
man scattering). Raman scattering can be classified 
into representations of the irreducible point group sym- 
metry of the crystal, selected by orienting the incident 
and scattered polarization vectors. The light scatter- 
ing polarizations allow energy fluctuations to be pro- 
jected onto different regions of the BZ, which has been 
useful in elucidating the behavior of electron dynamics 
for different momentum states [3. For example, for 
nonresonant scattering in tetragonal systems, the op- 
erator Als (q) = Ek £ ( k ) c k+ q /2 c k-q/2 = T(q) is the 
kinetic energy operator for q = Ai g scattering, and 

Csi g (q) = Ek e ( k + Q)4+ q /2 c k-q/2 is a modified ki- 
netic energy operator for q = Bi g scattering for the 
zone boundary wavevector Q = (tt, 0) [which is general- 
ized to higher dimensions as Q = (w, 0, n, 0, ...)]. Thus 
the Ai g scattering response is associated with energy 
fluctuations from all regions of the BZ, while the Bi g 
response is associated with dynamics near the BZ edges 
and away from the diagonals. The sum rule in Eq. J2J for 
nonresonant light scattering is similar to the optical sum 
rule, except the (model-dependent) operator average is 



different, and now depends on the potential energy since 
Tt = T+V [see Eq. (J2J] . In the general case, there are two 
contributions — one from the kinetic energy commutator 
and one from the potential energy commutator: 



(Wlmxm(q, v) = 7T^ Tr ( e 



-^[T(q,Q) + <7(q,Q)]). 

(5) 

The first term on the right-hand side of Eq. (JjjJ 
does not depend on the type of Hamiltonian and 
is given by f (q, Q) = - £ k (7 i(k, q, Q)c[ CT c k , CT , 
with i(k,q,Q) = e(k - Q - q/2) 2 '[e(k - q) - e(k)] + 
e(k- Q + q/2) 2 [e(k + q)-e(k)]. Regardless of the sym- 
metry [Q = for Ai g or (71", 0) for Bi g ] the kinetic- 
energy part of the sum rule does not contribute to the 
(q = 0) Raman scattering response, but it does con- 
tribute to the (q ^ 0) inelastic X-ray scattering response 
in a polarization-dependent manner. 

However the second term in Eq. 10 depends on the 
particular choice for the potential energy in the Hamil- 
tonian. Evaluating these sum rules for different mod- 
els is complicated, because the operator averages involve 
complex correlation functions. We consider the simplest 
model here, the spinless Falicov-Kimball model (FK) 



Hfk 



2Vd 



: c\c 3 + u 



(6) 



Here cj (cj) creates (destroys) a conduction electron at 
site i, uii is a classical variable (representing the localized 
electron number at site i) that equals or 1, t* is a renor- 
malized hopping matrix that is nonzero between nearest 
neighbors on a hypercubic lattice in (i-dimensions , U 
is the local screened Coulomb interaction between con- 
duction and localized electrons, and denotes a sum 
over sites i and nearest neighbors j. In our calculations, 
the average filling for conduction and localized electrons 
is set to 1/2. The FK model can be thought of as a limit 
of the Hubbard model for which one species of spin is not 
allowed to hop. 

For the FK model, the operator U in the sum rule is 



tW(q,Q) = ^ £ 4c i+s+5 ,B-^< s ^ 



2dN 



(7) 



i,6,S' 



x[e- l ?< s+s '\w i+s , - w+s+6>) ~ e^< s+s '\ Wi - w i+s )]. 

This operator involves the difference of correlated hop- 
ping operators, where the hopping of the conduction elec- 
trons is correlated with the presence of the localized elec- 
trons; note how the potential energy directly enters when 
5+8' = 0. The operator for the Hubbard model is similar, 
but involves complex spin and spin-flip hopping correla- 
tion functions, which could be measured, in principle, 
with neutron scattering. It is surprising that the sum 
rule for a charge response function can be related to a 
spin response function in the Hubbard model. 
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The result in Eq. J5J) is valid for any dimension. It 
can be evaluated exactly by examining the large dimen- 

I 



sional limit usin 
ory (DMFT) 



results from dynamical mean held the- 
Thc result for the FK model is 



dvv Im Xm (q, v) = ^Tr(T(q, Q) + [/(q, Q)> 



Im / dwf(w) \ (1 - X)X 



1 - -Z(u)G(u) - Z'\uj) [1 - Z{u)G(u)] 
f 1 — X' 2 

+Im / du>f(u)(E[u>) - U(wi)) I G{uj) + {Z{uj) - 1/G(lo)}X' 2 



+ -(A--l)[l-Z( W )G(w)] 



(8) 



with f(uj) = 1/[1 + exp(j3u>)] the Fermi function, G(uj) 
and the Green's function and self energy, respec- 

tively, and Z(ui) = u> + [i — The difference between 

A\ g and B\ g scattering comes from the momentum fac- 
tors X' = limrf^oo(l/d) 2i=i cos(qi/2)h l , with h = 1 
(— 1) for Ai g (Big) scattering. Note that both the hrst 
term associated with the kinetic energy and the second 
term associated with the potential energy contribute to 
the sum rule for inelastic x-ray scattering (q ^ 0) . On the 
other hand, only the second "potential energy" term con- 
tributes to the sum rule for Raman scattering (q = 0). 
Therefore, the sum rule's momentum dependence con- 
tains information regarding the potential and kinetic en- 
ergy contributions. 

Note that the sum rule for (X = 1) Raman scattering 
is proportional to the correlated part of the self energy, 
due to the factor U(wi) which has the Hartree shift 

subtracted off. The term in wavy brackets can be viewed 
as a Green's- function weighted variance of the bandstruc- 
ture for A\ g scattering, which can be interpreted as a 
measure of the width of the strongly-correlated band. 
Further, we note that in DMFT, where the self energy is 
local (and d — > oo), the sum rule for B\ g Raman scatter- 
ing would also apply to the second moment of the optical 
conductivity: 



/•oo />oo 

/ dvvImxB lg (q = 0, v) oc / dvv 1 a 
Jo Jo 



(v). (9) 



However, this relationship is violated if there is any mo- 
mentum dependence associated with the self-energy or in 
hnite dimensions. 

We plot the sum rule for the FK model in Fig. 2] for 
Bi g Raman scattering (q = 0) and different values of U 
as a function of temperature (results for A\ g are simi- 
lar). We have checked that an explicit integral of the 
nonresonant Raman response agrees with the sum rule 
plotted in Fig. ^ f° r many different values of U and 
T. We expect similar results to hold for other models 
of correlated electrons like the Hubbard model. Fig. 
shows that these sum rules are essentially constant at low 
temperature, indicating (i) the sum rule can be used to 




FIG. 1: Big sum rules as a function of temperature for the 
spinless FK model solved within DMFT (the Ai s results are 
similar). The different curves correspond to different values 
of the interaction strength increasing in steps of 0.1 from 0.5 
on the lowest curve to 4.0 on the highest curve. 



calibrate data from different samples and different tem- 
peratures; (ii) the sum rule can be used to determine 
the frequency above which interband transitions become 
prominent; (iii) the Raman response function multiplied 
by the frequency should be used to track spectral weight 
shifts because of the sum rule; and (iv) the sum rules 
have a momentum dependence for inelastic x-ray scatter- 
ing that can be calculated and compared to experiment. 
Note that unlike the f-sum rule, which is evenly weighted 
throughout the spectral range, the Raman sum rules are 
heavily weighted at higher energies due to the multipli- 
cation by the frequency. Consequently, unless there is a 
clear separation between low and high energy bands in 
the system, higher energy bands can significantly distort 
the sum rule. 

The momentum dependence of the sum rule is reduced 
as U increases and the physical properties become more 
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FIG. 2: Big scaled sum rules showing the scaling behavior 
with U 2 for the FK model discussed in the text. Deviations 
in scaling occur when U > U c — v2t*, the critical value of U 
for the metal-insulator transition in the half-filled FK model. 



local. The sum rule generally increases for increasing mo- 
mentum transfer since the kinetic energy term in Eq. © 
contributes to the sum rule as phase space is created for 
light scattering by increasing q. 

The underlying MIT is revealed via a simple scaling 
analysis, shown in Fig. El It can be shown that for metal- 
lic systems the sum rule varies as U 2 for small U , with 
large deviations occuring as the MIT is approached. In 
Fig. we plot the sum rule divided by U 2 for different 
temperatures as a function of U. We see that the sum 
rule data collapses for small U onto a single line for all 
temperatures. The scaled data abruptly fans out from 
the straight line at approximately U = U c and a strong 
temperature dependence emerges for larger U. Thus the 
onset of a strong temperature dependence and deviations 
from U 2 scaling can be used as a straightforward and 
quantitative way to identify a MIT from the optical data. 

We show the relevance of these sum rules to an experi- 
mental system that is ideal for this situation, SmB 6 . The 
Raman results on SmBg provide an ideal comparison to 
this model for several reasons: the spectra have a clear 
separation between the low energy (intraband) and high 
energy (interband) scattering contributions; only the low 
energy component in SmBg exhibits significant spectral 
weight changes due to the gap formation at low temper- 
atures; and the photon energy used in the experiment 
lies in a gap in the density of states of SmBg, and there- 
fore the Raman response is not likely to be influenced by 
resonant or mixed-scattering effects. When we take the 
Raman data, and multiply by the frequency, we find that 
the sum rule holds to within five percent in this system 
(see Fig.|3Jl, confirming the behavior determined from the 
FK model in DMFT (Fig. [TJ. 

In summary, we have discovered a class of sum rules for 
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FIG. 3: Raman response multiplied by frequency for SmE>6 at 
four different temperatures. On top is a plot of the integrated 
sum rule as a function of the frequency cutoff. Both curves are 
plotted in arbitrary units. The near crossing of the curves in 
the upper panel indicates where a cutoff can be placed (since 
the higher energy bands have little T dependence, the curves 
become parallel at high frequency, because the response is the 
same for all T). 



inelastic light scattering that are useful in analysing both 
Raman scattering and inelastic x-ray scattering. The 
sum rules depend crucially on the form of the electron 
interactions, and thereby yield useful information about 
the correlations in the material. While light scattering 
data on correlated metals or insulators is still rather lim- 
ited (particularly with regard to inelastic X-ray scatter- 
ing), our class of sum rules may be employed to analyze 
data and elucidate how the kinetic and potential energy 
of the system evolves as a function of doping, pressure, 
or temperature. 
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